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ABSTRACT 


t 

Strong  turbulence  may  be  produced  when  an  internal  gravity  wave 
encounters  a  critical  level.  This  turbulence  can  dissipate  the  wave  motion, 
so  that  the  momentum  transported  by  the  wave  is  deposited  into  the  mean 
flow.  Thus,  gravity  wave-critical  level  encounters  are  a  significant 
mechanism  by  which  the  wind  may  be  modified  by  gravity  waves. 

Such  an  encounter  has  been  computed.  The  rms  turbulent  velocities  are 
seen  to  be  a  substantial  fraction  (some  25%)  of  the  wave-induced  velocity. 
The  absorption  of  the  wave  is  seen  to  alter  the  wind  profile,  with 
turbulence  distributing  the  wave  momentum  flux  above  and  below  the  critical 
level.  A  variety  of  features  of  the  wave  and  turbulence  motion  may  be 
extracted.  One  which  will  be  useful  for  comparison  with  radar  observations 
is  the  refractive  index  structure  function  C£,  and  contour  plots  are 
displaced. 

Several  other  quantities  depicting  the  evolution  of  the  encounter  also 
are  presented.  Plans  for  a  detailed  analysis  of  these  and  new  computations 
are  outlined.  Also  discussed  are  comparisons  of  these  calculations  with 
Doppler  radar  observations  and  laboratory  measurements . 


INTRODUCTION 


The  specific  purpose  of  this  work  is  to  study  numerically  the  turbulent 
gravity-wave  critical-level  encounter.  The  larger  context  of  the  study  is 
the  manner  in  which  gravity  waves  alter  the  mean  wind.  As  will  be  seen  in 
results  presented  below,  this  encounter  can  be  a  site  of  strong  turbulence 
production  by  wave  breaking  as  the  wave  approaches  the  critical  level  (where 
the  wind  speed  in  the  direction  of  wave  propagation  matches  the  wave  phase 
speed).  This  turbulence  dissipates  the  wave  motion,  so  that  momentum 
transferred  by  the  wave  is  deposited  into  the  mean  flow. 

A  second  important  point  of  the  work  is  the  investigation  of  the 
characteristics  of  the  generated  turbulence.  Many  features  of  the 
turbulence  field  emerge  from  the  calculations;  one  of  the  most  interesting 
is  the  computation  of  the  refractive  index  structure  function  c£. 

This  quantity  may  be  deduced  from  the  strength  of  radar  backscatter  and  so 
is  a  natural  and  useful  point  of  convergence  between  the  numerical 
calculations  and  radar  observations  of  the  atmosphere.  These  results  are, 
to  my  knowledge,  the  first  calculations  of  their  kind. 

The  gravity-wave  critical  level  encounter  produces  strong  dynamical 
forces  and  rapid  transfer  of  energy  among  the  wave  spectral  components,  a 
reflection  of  its  well-known  vigorous  nature.  Because  of  this  vigor,  the 
numerical  calculations  have  been  found  to  be  rather  difficult  when 
turbulence  is  also  allowed  to  develop.  However,  an  emphasis  on  very  stable 
marching  algorithms  has  enabled  a  calculation  to  be  carried  out.  These 
computations  are  the  particular  important  achievement  of  the  work  for  the 
past  year.  Now  a  detailed  investigation  of  the  encounter  may  be  made  using 
these  and  new  calculations. 

The  next  section  highlights  these  results,  and  a  concluding  section 
outlines  the  course  of  the  study  for  the  next  year.  An  appendix  displays 
the  mathematical  equations  which  are  solved  numerically. 

NUMERICAL  CALCULATIONS 

In  order  to  minimize  ambiguity  in  interpreting  the  results,  the 
computation  has  been  arranged  to  model  an  ideal  case;  however,  the  values 
assigned  to  scales  and  parameters  are  certainly  commonly  realized  in  the 
atmosphere,  e.g.,  a  wave  encounter  with  a  shear  layer  associated  with  the 
polar  jet  stream. 

The  wind  is  described  by  the  frequently  used  hyperbolic  tangent  profile 


u(z)  =  U  Itanh(z/H)  ♦  1}, 


and  the  potential  temperature  is  exponentially  increasing  with  height  z  so 
that  the  Brunt-Vaisala  frequency  is  constant: 


N(z)  *  N. 
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The  values  assigned  to  the  velocity  and  height  scales  are 


U  s  14.142  meters/sec 
H  =  1000  meters 


and 


N  »  0.0112  rad/sec. 


Hence  the  minimum  Richardson  number  is 


Rimin  *  (NH/U)2  *  2. 


A  plot  of  the  wind  profile  is  shown  in  Figure  1.  The  height  range  shown  is 
that  of  the  calculation: 


-  4  <  z/H  <  2. 


In  this  figure  and  all  following  figures,  all  quantities  are  non-dimension- 
alized  by  U  and  H.  The  exception  is  c£,  which  is  given  dimensionally 
(i.e.,  in  meters-2^). 

The  wave  is  forced  at  the  bottom  of  the  computational  domain  by  a 
boundary  condition  equivalent  to  the  requirement  on  the  wave-induced 
vertical  velocity  component  w  that 


w(x,z,t)  =  F(t)  W  sin(kx  -  ut),  z/H  =  -  4, 


where 


W/U  - 

0.1 

kH 

0.5 

wH/U  a 

0.5 

The  amplitude  function  is 

F(t)  =  sin2  (w  t/T) 
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with 


T  -  3  wave  periods. 


Hence  the  wave  emerges  with  initially  small  but  increasingly  larger 
amplitude  (avoiding  startup  transients),  and  after  three  periods  propagates 
away  from  the  bottom  boundary  without  being  re-generated  there.  About  one 
wave  period  is  generated  with  substantial  amplitude.  Although  the 
generation  amplitude  W  is  not  small,  neither  is  it  unusually  large;  that  is, 
the  dramatic  evolution  of  the  wave  as  it  approaches  the  critical  level  is 
not  due  to  a  synthetically  enormous  wave  amplitude,  but,  rather,  to  the 
fluid  dynamics  associated  with  such  encounters. 

Note  the  wave  parameters  are  chosen  so  that  the  critical  level  is  at 
z  =  0: 


u(z)/U  =  (wH/U) / (kH)  at  z  =  0. 


The  positive-definite  turbulence  correlations  are  set  to  the  value  10~7, 
which  essentially  is  zero.  The  correlations  (e.g.,  momentum  flux)  are  set 
exactly  to  zero. 

The  number  of  grid-points  in  the  z  direction  is  151,  and  there  are  24 
spectral  modes,  or,  with  the  Fast  Fourier  Transform  conventions,  48  points 
in  the  horizontal.  The  calculation  extends  to  4.25  periods,  with  a  time 
step  size  of  0.01  periods. 

An  overview  of  the  evolution  of  the  encounter  is  obtained  from  Figure  2. 
Plotted  there  are  the  minimum  value  of  the  total  (mean  plus  wave)  potential 
temperature  gradient,  maximum  turbulent  variance  in  the  horizontal  velocity 
component,  and  maximum  turbulent  heat  flux.  The  potential  temperature 
gradient  initially  decreases  from  2  due  to  emergence  of  the  wave  at  the 
lower  boundary;  the  wave  amplitude  there  is  not  sufficient  —  even  at  the 
maximum  at  t  =  1.5  periods  —  to  drive  this  gradient  negative.  However,  as 
the  wave  encounters  the  critical  level,  there  the  larger  Lagrangian 
displacements  of  fluid  results  in  denser  air  being  pulled  up  over  less 
dense;  this  is  indicated  by  the  temperature  gradient  becoming  negative. 

About  one-half  wave  period  later  the  instability  of  this  arrangement 
produces  a  turbulent  heat  flux  which  almost  immediately  leads  to  rapid 
generation  of  turbulent  velocity  fluctuations.  The  amplitude  of  the 
variance  shown  in  Figure  2  is  about  0.06,  corresponding  to  a  rms  value  of 
about  0.25.  Since  the  (nondimensional)  horizontal  wave-induced  velocity 
maximum  is  about  1.0,  we  see  that  the  turbulent  velocity  fluctuations  are 
quite  significant  compared  to  the  wave  velocities.  As  the  turbulence 
restores  the  stable  stratification,  its  source  of  energy  is  correspondingly 
diminished,  and  eventually  the  turbulent  motions  die  away  under  the  combined 
influence  of  dissipation,  diffusion,  and  stratification. 
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This  encounter  is  more  vividly  depicted  by  plots  in  the  x  -  z  (i.e.t 
computational)  plane.  Accordingly,  contour  plots  of  potential  temperature, 
turbulent  kinetic  energy,  and  C^,  and  vector  plots  of  the  wave-induced 
velocity  are  shown  at  time  1.5,  3.0,  3.5,  and  4.25  periods  in  Figures  3,  4, 
5,  and  6,  respectively.  In  these  plots,  the  z  coordinate  again  ranges  over 
-  4  <  z/H  <  2,  and  the  x  coordinate  ranges  over  on  wavelength  *  4*. 

In  Figure  3,  the  upper  edge  of  the  wave  is  just  reaching  the  critical 
level  at  z  *  0.  Most  of  the  wave  is  still  in  the  region  of  little  shear, 
the  wave  phase  lines  are  nearly  straight,  and  the  displacements  revealed  by 
the  shape  of  the  potential  temperature  contours  are  approximately  sinusoidal 
(linear).  No  turbulence  has  been  generated. 

In  Figure  4,  at  3.0  periods,  the  wave  is  well  into  its  encounter  with 
the  critical  level.  Just  below  this  level,  the  potential  temperature 
contours  display  the  pronounced  'overhang*  of  denser  air  over  less  dense 
air.  The  vector  plot  shows  that  the  wave  is  just  beginning  to  plunge 
downward,  driven  by  the  buoyancy  forces  associated  with  this  overhang.  This 
unstable  stratification  also  has  resulted  in  an  almost  explosive  growth  of 
turbulence;  note  the  coincidence  of  the  turbulent  kinetic  energy  with  the 
region  of  unstable  stratification.  There  has  been  some  discussion  in  the 
literature  about  whether  these  encounters  will  generate  turbulence  by  this 
sort  of  convective  instability,  or  if  the  mechanism  is  through  shear  (or 
'dynamic')  production.  This  result  unequivocally  indicates  that  it  is  the 
former  case  for  waves  of  this  period  and  wavelength.  The  contour  plot  of 
Ca  shows  that  it  generally  follows  the  gross  features  displayed  by 
the  contour  plot  of  turbulent  kinetic  energy,  as  may  be  expected.  The  value 
of  mean  density  used  in  computing  is  a  tropopause  value  of  0.42  kg/m. 

The  peak  of  in  this  plot  is  about  10-14  meters  ~a/*.  which  is 
roughly  two  orders  of  magnitude  greater  than  that  observed  by  Doppler  radar 
during  wave  events  at  this  altitude.  However,  the  geometrical  extent  of  the 
peak  values  is  rather  small,  much  less  than  a  radar  beamwidth,  so  the 
calculations  seem  satisfactory. 

By  3.5  periods  (Figure  5)  the  downward  plunge  of  the  wave  has  taken 
place  and  energy  has  cascaded  to  higher  spectral  modes,  as  seen  from  the 
vector  plot  given  in  this  figure.  These  higher  modes  are  evanescent,  not 
propagating  away.  Much  of  the  original,  identifiable  structure  of  the  wave 
has  been  destroyed.  The  relatively  strong  stable  stratification  so  created 
inhibits  the  downward  diffusion  of  the  turbulent  kinetic  energy  and  results 
in  a  concentration  of  the  refractive  index  fluctuations  (as  indicated  by  the 
contour  plot  of  C^)  near  the  bottom  of  the  region  of  turbulent  activity. 
These  features  are  qualitatively  similar  to  those  seen  in  photographs  of 
laboratory  gravity-wave  critical  level  encounters.  Remnants  of  the  wave 
motion  advect  the  turbulence  upward  into  levels  of  greater  wind  speed  (e.g. , 
near  z  *  0),  and  this  process  is  aided  by  upward  diffusion.  The  greater 
wind  then  advects  the  turbulence  horizontally,  stretching  it  out  nearly 
across  the  entire  wavelength. 

Figure  6  displays  the  situation  at  4.25  periods.  Little  of  the  form  of 
the  original  wave  motion  remains;  the  wave  motion  is  dominated  by  the 
'sloshing'  of  the  shorter  wavelength  evanescent  modes.  Notice  that  in  the 
region  of  remaining  turbulent  activity  there  is  little  wave  motion,  it 
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having  been  damped  out  by  the  turbulence.  The  turbulence,  under  the 
influence  of  stable  instead  of  unstable  stratification,  is  rapidly  dying 
away.  The  different  advection  rates  by  wind  at  differing  levels  has  formed 
the  remaining  turbulence  into  a  fairly  thin  quasi-homogeneous  layer. 

Finally,  the  wind  profile  at  this  last  time  of  4.25  periods  is  shown  in 
Figure  7.  Even  though  the  wave  was  generated  with  substantial  amplitude  for 
only  about  one  period,  the  significant  increase  in  wind  speed  around  the 
critical  level  is  evident. 

RESEARCH  PLANS  FOR  THE  COMING  YEAR 

In  the  context  of  modification  of  the  wind,  a  central  point  of  inquiry 
is  the  behavior  of  the  momentum  transport  by  the  wave  motion  and  the 
turbulence  as  the  critical  level  encounter  proceeds.  Accordingly,  a  focus 
of  the  upcoming  work  will  be  on  this  transport,  as  revealed  by  the 
calculations,  in  particular,  on  how  the  momentum  transport  is  shared  between 
the  wave  motion  and  the  turbulence,  and  how  it  is  distributed  over  the  mean 
flow  as  the  wave  is  absorbed.  Further,  the  heat  flux  plays  a  major  role  in 
the  critical  level  encounter,  as  may  be  seen  in  the  results  presented 
above.  Thus,  this  transport  also  will  be  an  important  subject  of  the 
study.  There  have  been  in  recent  years  several  speculations  in  the 
literature  on  the  relative  magnitude  of  the  momentum  and  heat  transports  in 
wave  breaking  events,  usually  cast  in  terms  of  the  question  as  to  the 
magnitude  of  the  Prandtl  number.  The  calculations  presented  above  and  those 
to  be  performed  allow  an  immediate  and  direct  computation  of  the  magnitudes 
in  question. 

These  mean  momentum  and  heat  fluxes  mediate  the  modification  of  the  mean 
wind  and  temperature  profiles  by  the  breaking  wave.  The  fluctuating 
components  of  these  .axes  modify  and  dissipate  the  wave  itself.  This 
direct  interaction  between  the  wave  motion  and  the  turbulence  is  an 
important  feature  in  its  own  right,  and  will  be  an  important  point  of 
investigation  in  the  upcoming  study. 

The  new  calculations  will  include  cases  where  the  wave  is  forced  for 
more  than  the  approximately  one  wave  period  of  the  case  shown  above.  Such 
longer  wave-train  cases  are  more  typical  of  atmospheric  waves. 

Plots  of  C*  such  as  those  given  above  will  be  made  for  all  cases 
computed.  Radar  records  in  which  suitably  monochromatic  wave  events  have 
been  identified  will  be  compared  to  these  calculated  results.  The  Poker 
Flat,  Alaska  and  the  Sunset  Colorado  Doppler  radars  appear  to  be  the  best 
candidate  instruments  for  such  comparison.  Even  order-of -magnitude 
agreement  will  lend  confidence  in  the  atmospheric  relevance  of  these 
calculations . 

Finally,  the  six-foot  diameter  annular  tank  constructed  and  operated  at 
Northwest  Research  Associates  (formerly  Physical  Dynamics,  Inc.)  may  provide 
measurements  of  breaking  waves  which  will  enable  fairly  detailed  comparison 
with  results  of  these  numerical  calculations.  Such  comparison  would  both 
test  the  numerics  as  well  as  provide  insight  into  the  breaking  wave 
observations . 
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APPENDIX:  MATHEMATICAL  FORMULATION 


a  <u;  u;> 
at~ 


j  (<u;  u;>,  <v>)  +  (cq  -7-  <u'k  u;>  a,  <u;  u;>) 

2  , 

+  P,,  +  G„  +  i ,  “  3  E 


q  2  y  <u;u;> 

p()  2  -<u;  u;>  ak  <u,>  -  <u;  u;>  ak  <u,> 

G,,  s  p,  <0'  u',>  +  0,  <9'  u',> 

i,  s  in  +  *..2  +  i,3 

2 

■  -C,  y  (<u;u;>  -  3-  d„q) 

f,2  s  -c2(P„  -  T^P).ps  T  P- 

♦,|3  =  ~c3  (G.,  -  -3  d,,  G)  ,  G  2  y  G„ 
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Figure  2.  EVOLUTION  OF  MINIMUM  POTENTIAL  TEMPERATURE  GRADIENT,  MAXIMUM 
HEAT  FLUX,  AND  MAXIMUM  VARIANCE  IN  HORIZONTAL  VELOCITY. 
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Figure  AT  t  «  3.0  PERIODS,  CONTOUR  PLOTS  OF  POTENTIAL  TEMPERATURE,  TURBULENT 
KINETIC  ENERGY,  AND  C£,  AND  A  VELOCITY  VECTOR  PLOT. 
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